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■ 

■ We study families of one-dimensional matter-wave bright solitons supported by the competition 
' of contact and dipole-dipole (DD) interactions of opposite signs. Soliton families are found, and 
, their stability is investigated in the free space, and in the presence of an optical lattice (OL). Free- 
• space solitons may exist with an arbitrarily weak local attraction if the strength of the DD repulsion 

■ ' is fixed. In the case of the DD attraction, solitons do not exist beyond a maximum value of the 

local-repulsion strength. In the system which includes the OL, a stability region for subfundamental 
, solitons (SFSs) is found in the second finite bandgap. For the existence of gap solitons (GSs) under 

( ' the attractive DD interaction, the contact repulsion must be strong enough. In the opposite case of 

the DD repulsion, GSs exist if the contact attraction is not too strong. Collisions between solitons 

1 I , in the free space are studied too. In the case of the local attraction, they merge or pass through 

each other at small and large velocities, respectively. In the presence of the local repulsion, slowly 
Qh ' moving solitons bounce from each other. 

d ■ 

.^H . PACS numbers: 03.75.Lm; 32.10.Dk; 05.45.Yv 

'■ I. INTRODUCTION 

> ■ 

^■f^ ' Stable localized matter-wave structures in Bose-Einstein condensates (BECs) are supported by the inter- 

fT^ , play between the intrinsic nonlinearity, which is induced by collisions between atoms, quantum pressure, 

OO ' which originates from the kinetic energy of atoms, and external potentials This mechanism has made 

I it possible to create bright solitons in condensates of ^Li and ^^Rb atoms confined in cigar-shaped traps 

T^lj- i 0, where the inter-atomic interactions are made attractive by means of the Feshbach resonance (FR) [3]. 

' In the condensate of *^Rb atoms with repulsive interactions, the introduction of an optical-lattice (OL) 

0^ . potential gives rise to gap solitons (GSs), as demonstrated experimentally in Ref. Q, see also review [H]- 

' Dark solitons have been created too, by means of various techniques, in the self-repulsive ^"^Rb condensate 

^ , 0. In terms of the theoretical description, the limit case of a very deep OL may be mapped, by means of 

r ■ the tightly-binding approximation, into a discrete nonlinear Schrodinger equation and, accordingly, GSs are 

r> I mapped into staggered discrete solitons 

■ New possibilities for the formation of matter-wave solitons are suggested by the presence of long-range 
" " ' interactions in dipolar condensates, which may be composed of magnetically polarized ^^Cr atoms @, dipolar 

molecules [9], or atoms in which electric moments are induced by a strong external field [lo| . Solitons 
supported by the dipole-dipole (DD) interactions were predicted in two-dimensional (2D) settings. In the 
isotropic configuration, with moments fixed perpendicular to the plane, the natural DD interaction gives rise 
to repulsion, which can support delocalized states in the form of vortex lattices In principle, the 

sign of the DD interaction in this configuration may be reversed by means of rapid rotation of the dipoles 
[l3 | , suggesting a possibility to create isotropic solitons [3| , as well as solitary vortices jl^, |15|] . On the 
other hand, stable anisotropic solitons have been predicted assuming the natural DD interaction between 
dipoles with a fixed in-plane polarization [l^. In addition to these results pertaining to the BEC context, it 
is relevant to mention that stable vortex rings were predicted in an optical model with the nonlocal thermal 
nonlinearity [iTj , and , elliptically shaped spatial solitons were created in such media experimentally [18] . 

Although one-dimensional (ID) configurations may be simpler than their 2D counterparts, ID matter- 
wave bright solitons were not yet studied in detail in models of dipolar condensates, except for discrete 
solitons of the unstaggered type, which were recently predicted in the condensate trapped in a deep OL 
[l^ . In those works, both attractive and repulsive signs of the onsite (contact) nonlinearity and long-range 
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DD interactions between sites of the respective lattice were considered. The objective of the present work 
is the theoretical study of various types of bright solitons possible in the continuum BEC model featuring 
the competition between the contact and DD interactions. The effective strength of the DD interactions in 
the ID geometry can be controlled by adjusting the orientation of the dipoles with respect to the axis of 
the linear trap, while the strength of the contact interactions may be effectively tuned by means of the FR 
technique, as shown in the condensate of ^^Cr atoms [lOl- The model is formulated both in the free space 
(i.e., in the absence of external potentials) and in the presence of the OL potential, which opens additional 
possibilities for the creation of stable localized states, including GSs (note that staggered discrete solitons, 
which, as mentioned above, correspond to GSs in the deep-OL limit, were not considered in Ref. [l^). In 
the framework of the Gross-Pitaevskii equation (GPE) with the ordinary local nonlinear term, the concept 
of GSs was elaborated in detail, see Refs. [2l[-[23| and review f^. However, to the best of our knowledge, 
it was not yet extended to the new physically relevant case, when the OL potential acts together with the 
long-range DD interaction - a situation that we address below. 

The paper is organized as follows. The model is formulated in Section II, which, in Section III, is followed 
by the consideration of solitons supported by competing nonlinearities - attractive local/repulsive DD, or 
vice versa - in the free space. The model which combines the competing nonlinearities of both types and 
the OL is considered in Section IV. In that case, we report results for regular solitons in the semi-infinite 
gap, and for GSs in the two lowest finite bandgaps. In Section V, we deal with collisions between moving 
solitons in the absence of the OL. The paper is concluded by Section VI. 



II. THE MODEL 



Our aim is to construct soliton states within the framework of the ID GPE for the mean-field wave 
function, ■0(a:,t). The scaled equation includes the OL potential, V{x) = esin^a; (where e is the strength of 
the OL, while its period is normalized to be tt), the local nonlinear term with the respective coefficient, g,,, 
and its nonlocal counterpart, with coefRcient gd, which accounts for the DD interactions: 
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where subscripts denote partial derivatives. The kernel K of the DD term was taken in two different ways, 
so as to avoid the divergence at x = x' . The first version relies on an explicit cut-off {CO kernel), 
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where yc is a constant. The other choice makes use of the regularized expression deduced in Ref. by 
means of the single-mode approximation (SMA), i.e., the SMA kernel: 
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with evic{y) being the standard complimentary error function. The two kernels are compared in Fig. [TJ The 
most significant difference between them is that the SMA version features a cusp at y = 0, whereas the CO 
kernel has a smooth maximum. We have concluded that the choice of yc — in CO expression 

which makes areas beneath both curves equal, provides for the best proximity of the corresponding results 
to what has been found using the SMA approximation. Below, we report results obtained with the SMA 
kernel, as its CO counterpart yields virtually identical findings. 

We fix the normalizations in Eq. ([T]) by setting gd = ±1, and then vary gc- In the the case of the ^^Cr 
condensate, a characteristic value of the relative strength of the DD and contact interactions, which may be 
estimated as \gd/9c\: is — 0.15 [l^l; this value may be altered in broad limits by means of the FR technique. 
The interactions are repulsive or attractive for gd,gc > and gd,gc < 0, respectively. We will focus on 
the case of competing interactions, with gcgd < 0, which is the most interesting one; in the case when both 
nonlinearities have the same sign, results turn out to be very similar to those reported previously in the 
local model. It is relevant to mention that, in terms of discrete systems, the competition of on-site (local) 
and inter-site (short-range nonlocal) interactions in ID and 2D Salerno lattices were considered in Refs. [2^ 
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FIG. 1: (Color online) The solid and dashed lines show the cut-off (CO) kernel, with yc = n and its single-mode 
approximation (SMA) counterpart, which are based on Eqs. ([2)) and ([3]), respectively. 



and [2a], respectively. A number of stable discrete-soliton states which are impossible in the absence of the 
competition were reported in those works, including cuspons and peakons in ID, and vortex breathers in 
2D. 

Stationary solutions to Eq. ([T]) with chemical potential /i are sought as ili{x,t) — 'I'(x) exp(— i/it). To 
construct such solutions, we discretize the resulting equation for "^{x) by means of a finite-difference scheme, 
which leads to a set of coupled algebraic equations. 
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with Xn = nAx, and Kn-m = K{\n — m|Aa;). Solutions to Eq. ([4]) are sought by dint of the Newton- 
Raphson scheme. Results were obtained with a reasonable accuracy by choosing Ax — 7r/40. To present 
soliton families, we will use norm N and width W of the soliton. 
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These definitions were adapted to the finite-difference form of the model as well. 

Stability of the solutions was analyzed in the standard way (see, e.g., Ref. ^^), by considering a pertur- 
bation in the form of 5ip{x,t) = exp(— i/it)[P(x) exp(— iAi) + Q* (x) exp{iX*t)] (where * stands for complex 
conjugate). The linearized equations for the perturbation eigenmodes (i.e., the Bogoliubov - de Gennes 
equations) are written as 
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The instability sets in when there emerges an eigenvalue with Im(A) ^ 0. Stability eigenvalues were obtained 
from a numerical solution of Eq. ([6]), by means of a standard eigenvalue solver from Fortran-based software 
package LAPACK [28.]. 



III. SOLITONS IN THE ABSENCE OF THE OPTICAL LATTICE 



We start the presentation of results by considering the model with competing nonlinearities in the free 
space, i.e., with e = in Eq. ([T]). This implies that solitons may exist with /i < 0. 



A. Attractive local and repulsive nonlocal interactions 



In accordance with what was said above, we first fix = 1 (the repulsive DD interaction), and vary 
negative Qc (local attraction) and negative /z. In this case, solitons exist for every Qc < 0; however, due 
to discretization, the numerical solution of Eq. ([l]), with the above-mentioned choice of Aa; = 7r/40, yields 
solitons in the region of \gc\ > 0.25. In fact, although very narrow solitons indeed exist at arbitrarily small 
values of —gc, the soliton's width becomes comparable to or smaller than this value of Ax at \gc\ ^ 0.9, 
i.e., a better numerical accuracy is required to produce accurate soliton solutions in this range. Note that, 
for very narrow solitons with amplitude 5*0, the nonlinear part of Eq. with g^ = 1, takes the form of 
{gc + KqAx) ^fp. The existence of the soliton demands a negative coefficient in this expression, i.e., 

\gc\ > (lOA)Ax, (10) 

where it was taken into regard that Kq = IO/tt, as per Eq. ([3]). In particular, for Ax = 7r/40, condition (fTO|) 
amounts to j^d > 0.25, as said above. 

Figure dja) shows typical examples of the solitons in the present case, while Fig. [3] represents soliton 
families in terms of dependences N{^) and W {gc), for fixed values of gc and /i, respectively; note that plots 
in panel (b) of the latter figure are cut at gc = —1, as the numerical accuracy is insufficient to extend 
them to smaller values of \gc\, as explained above. The stability of the solitons was verified both through 
the computation of the eigenvalues, using Eq. and by means of direct simulations of the evolution of 
perturbed solitons. The well-known Vakhitov-Kolokolov criterion, dN/djjL < |29|], also suggests the stability 
of solitons in this case, although the negative slope of the N{jjL) curve in Fig. [Slja) is very small. 



B. Repulsive local and attractive nonlocal interactions 



Now, we fix = — 1, varying > and /i < 0. Contrary to the previous case, solitons (which are stable) 
can be readily found for all values of gc up to gc ~ 3.6, the difference between gc = (zero local interaction, 
while the DD attraction is present) and gc > amounting to a gradual increase of the soliton's amplitude 
and width with gc- As seen in Fig. [U^b), close to gc — 3.6 the soliton develops a compacton-like shape, and 
solitons cannot be found at gc > (5c)i„ax ~ 

The existence of (gc)inax '-^^^ easily explained. Indeed, in the limit of Ax — > 0, and for a very broad 
soliton, the nonlinear part of Eq. (U), with gd = —1, takes the approximate form of 
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The necessary condition for the existence of solitons is that the coefficient in front of \tlj{x)\'^'ip{x) in this 
expression must be negative [cf. the derivation of Eq. pop ], i.e., 

/ + 00 
K{y)dy = 20^-3/2 « 3.59, (12) 
-OO 



5 




-3 -2 -1 1 2 3 -40 -20 20 40 

X X 



(a) (b) 

FIG. 2: (Color online) (a) Profiles of stable solitons in the model with attractive local and repulsive nonlocal inter- 
actions. The solitons were found in the numerical form with resolution Ax = 7r/40, which yields accurate results for 
\gc\ > 0.9. The peak gets sharper as [gd decreases, (b) Profiles of solitons in the model with repulsive local and 
attractive nonlocal interactions, obtained with Ax = tt/IO. The compacton-like solution corresponding to \gc\ = 3.6 
is close to the critical point, beyond which (at larger \gc\) solitons do not exist, irrespective of Aa;. In both panels, 
all solitons pertain to fj, = —1. 



where expression ([3]) was used to perform the integration. For finite Ax, the integral in expression (jl2p is 
replaced by X^m -^m^^;. In particular, for Ax = tt/IO, this yields (5c),„ax ~ 3.708, which is consistent with 
the above-mentioned numerical finding. Note also that the vanishing of the coefficient in front of \ip{x)\'^'ipi^) 
in expression (jlip at (fl'c)iiiax ~ 5c — > implies that the soliton's amplitude and norm diverge in this limit, 
which is corroborated by the numerical results shown in Fig. 21 

The dependence of the soliton's norm and width on the chemical potential is displayed in Fig. [5l Unlike 
the nearly fiat N{fi) dependences in Fig. [31 the present ones clearly satisfy the VK stability criterion, 
dN/d^ < 0. The full stability of the solitons was confirmed by the computation of eigenvalues, using Eq. 
®. 

At this point, it is also worth to briefiy consider dark solitons, which are known to be stable in BEC with 
contact repulsive interactions, such as *^Rb condensates 0,0]. Dark solitons in three-dimensional dipolar 
BECs were recently considered in Ref. [lOl, where it was shown that, for sufficiently strong repulsive DD 
interactions and a sufficiently deep OL in the soliton's nodal plane, dark solitons exist and are stable. In the 
present ID setup, to investigate the existence and stability of dark solitons, it is first necessary to ensure that 
the respective background, namely, constant-amplitude state 4' = y/l^Gxp{—i^t), is modulationally stable. 
A comprehensive analysis of the modulational instability (MI) of the background in the context of Eq. ([T]) 
can be performed following the lines of Ref. (3]| . Here we will briefiy consider this issue and provide an 
example of a stable dark soliton, assuming gc + {20/TT^^^)gd > 0. In this case, the effective nonlinearity for 
long- wavelength perturbations (i.e., those with wave numbers /c ^ 0) is self-defocusing, i.e., the MI band 
cannot start at fc = 0, as it does in the case of the standard nonlinear Schrodinger equation with the self- 
focusing nonlinearity. Although the MI band may appear at finite k, this is not expected to happen as long 
as the background density is small enough, because the maximum MI gain is proportional to that density. 
Thus, in this case, a modulationally stable background may exist, and dark-soliton solutions can be found. 
As an example, in Fig. |S1 we show a stable dark soliton (its stability was verified through the computation 
of the full spectrum of eigenvalues for small perturbations), which was found for = 1 and gc = 5, gd = —I- 
A systematic analysis of the MI of the background and dark-soliton families in the framework of Eq. ([T]) is 
beyond the scope of this work, and is deferred to a separate publication. 



6 



-1.8 -1.6 -1.4 -1.2 



-2.5 



-0.8 -0.6 -0.4 -0.2 



0.4 r 



0.35- 



0.3- 



0.25- 



0.2- 



0.15- 



0.1 



H=-2 



-2.5 -2 



-1.5 



(a) 



(b) 



FIG. 3: (Color online) (a) The norm of solitons in the free space, in the case of attractive local and repulsive nonlocal 
interactions, versus the chemical potential, (b) The width of the solitons, defined as per Eq. ([5}, versus the strength 
of the local attraction. In panel (b), the plots are shown for W > Ax = 7r/40: recall that obtaining numerical results 
for smaller W requires using smaller Aa;. 
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FIG. 4: (Color online) The norm and width of fundamental solitons in the free space, in the case of attractive nonlocal 
and repulsive local interactions, versus the strength of the local repulsion. 



IV. SOLITONS IN THE OPTICAL LATTICE 



In the presence of the OL potential, generic results for regular solitons and gap solitons (GSs) in the model 
with the competing interactions can be adequately represented by fixing the OL strength to e = 6, which 
is adopted below. GSs have been found in the first and second finite bandgaps of the OL-induced linear 
spectrum. For e = 6, the two numerically computed (with Ax = 7r/40) bandgaps cover, respectively, the 
following intervals of the chemical potential: 



1.61 < iK 4.26 and 4.63 < ^ < 6.02. 



(13) 
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FIG. 5: (Color online) The norm and width of fundamental solitons in the free space, in the case of attractive nonlocal 
and repulsive local interactions, versus the chemical potential. 
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FIG. 6: (Color online) The profile of a stable dark soliton in the model with attractive local and repulsive nonlocal 
interactions, for gc = 5, gd — —1, and /i = 1. The soliton was found in the numerical form with resolution Ax = tt/IO. 



A. Attractive local and repulsive nonlocal interactions 

Here, we consider solitons in the case of the competition between local attraction {gc < 0) and repulsive 
DD interactions {gd = 1), varying /i and gc- The self-focusing character of the local interaction allows the 
existence of regular solitons in the semi-infinite gap, which is —oo < fi < 1.5810 for e — 6. The numerical 
solution of Eq. with Ax = 7r/40, yields regular solitons for gc < —0.25, similar to the case of e = 0, see 
above. 

Apart from the solitons in the semi-infinite gap, GSs have been found in parts of the first and second 
finite bandgaps, for gc exceeding a certain critical value, which depends on /i, as shown in Fig. [T] Multi- 
humped solitons, which are bound states of fundamental single-humped solitons, can be found too. Numerical 
results demonstrate that the existence range for the multi-humped solitons is slightly broader than for the 
fundamental ones. In the second bandgap, single-humped solitons exist (and are stable) at /i > /id- ~ 5.1. 
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FIG. 7: (Color online) Unshaded and light-gray-shaded areas represent existence regions for solitons in the semi- 
infinite and two lowest finite gaps, in the model combining the optical lattice with competing attractive local and 
repulsive nonlocal interactions. Dark-gray shading covers Bloch bands, where solitons do not exist. In the gray- 
shaded parts of the finite bandgaps (intermediate between light and dark gray), the solitons do not exist either. In 
the white part of the second bandgap on the left side of the vertical dashed line, only subfundamental solitons (SFSs) 
exist. They are stable between the red (solid) line and the gray area. On the right side of the dashed line, stable 
fundamental solitons exist in the white region, while unstable SFSs exist both in the white and light-gray areas. 

This boundary value, which corresponds to the vertical dashed hnc in Fig. [71 is located above the lower edge 
of the second bandgap, fi « 4.63, see Eq. (fTS]). 

Another species of stable GSs was found in a part of the second finite bandgap, in the form of subfunda- 
mental solitons (SFSs). These are antisymmetric modes which are squeezed, essentially, into a single cell of 
the OL. The norm of the SFS is smaller than that of the fundamental GS, if the latter one can be found at 
the same value of /x (hence the name of "subfundamental" [32l] ) . In a narrow interval adjacent to the lower 
edge of the second bandgap, 4.63 < fj. < 4.71, the SFSs are stable (see Fig. [Tj), whereas above this interval, 
they undergo a Hamiltonian-Hopf bifurcation, being unstable for \gc\ smaller than a certain critical value. 
In the region of /i > ficr ~ 5.1, SFSs are unstable in their entire existence region. In direct simulations, the 
destabilized SFSs spontaneously transform themselves into stable fundamental solitons belonging to the first 
(rather than second) bandgap; a similar scenario of the instability development of SFSs was found in the 
local model [111 . 

The existence regions for fundamental solitons (both regular ones in the semi-infinite gap and GSs in the 
first two finite bandgaps) and SFSs in the {fJ.,gc) plane are depicted together in Fig. [71 the existence region 
for multi-humped bound states is not shown separately, as it almost coincides with that of the fundamental 
solitons. In addition. Fig. [5] shows the norm as a function of /i at f/c = — 1, for families of subfundamental 
and fundamental GSs, and for stable bound states of fundamental GSs [in the semi-infinite gap, the norm 
of all types of regular solitons very weakly depends on /i, cf. Fig. [H^a)]. Typical examples of solitons of all 
these types are displayed in Fig. [51 Bound states in the second finite bandgap are not shown, as they are 
completely unstable against oscillatory instabilities (while they are stable in the semi-infinite and first finite 
gaps). In fact, the same instability of bound states of GSs in the second finite bandgap occurs in the local 
model. 



B. Repulsive local and attractive nonlocal interactions 

To adequately represent results in the model featuring the competition between the local repulsion (gc > 0) 
and DD attraction (g^ = —1) in the model with the OL, it was sufficient to use a coarser numerical mesh, 
with Aa: = n/lO (recall that Ax — 7r/40 was used above). The nonlocal self-attraction allows the existence 
of solitons in the semi-infinite gap. Similar to what was reported above for the same case in the free-space 
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FIG. 8: (Color online) N{fi) curves for gap solitons in the model including the optical lattice, local attraction, and 
nonlocal repulsion. The strength of the local attraction is fixed to Qc — —1. 
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FIG. 9: (Color online) Profiles of stable solitons labeled A-F in Fig. |8l Additionally, panels G-J display examples 
of stable regular solitons found in the semi-infinite gap. The parameters are (qc — —1.5, jj, — 1) and (qc = — 1, 
jj, — 3 and 5.5) for the solitons in the semi-infinite gap, and first and second finite bandgaps, respectively. 



model, the solitons become very broad at gc approaching the critical value, gc ~ 3.6. Other features of the 
regular solitons found in the semi-infinite gap are also similar to those of their free-space counterparts. The 
similarity holds also in the case of gc = 0, i.e., in the model with the pure nonlocal attractive interactions. 

The existence range for stable fundamental solitons (both the regular ones and GSs) and SFSs in the 
present version of the model is depicted in Fig. [TOl Further, Fig. [TT] shows the respective N{^) dependences, 
including those for bound-state solutions. This figure shows iV(/i) lines in the semi-infinite gap too, as, on 
the contrary to the situation for the model with gc < and gd = 1, these lines do not degenerate into 
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FIG. 11: (Color online) A^(/i) curves for regular (left panel) and gap (right panel) solitons in the model with the 
optical lattice, local repulsion, and nonlocal attraction. The strength of the local repulsion is fixed to Qc ~ 1 (a) and 
5c = 5 (b). 



N — const. In fact, the opposite signs of dN/dfi for regular solitons and GSs is a typical feature, observed 
in local models as well. Typical profiles of various soliton species belonging to the semi-infinite and finite 
gaps are displayed in Fig. [121 



V. COLLISIONS BETWEEN MOVING SOLITONS 



We have also studied collisions between solitons moving in the free space. In the model with the repulsive 
local and attractive DD interactions, a usual collision scenario is observed: at small velocities, solitons merge 
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FIG. 12: (Color online) Soliton profiles at points A-J in Fig. 1111 The parameters are (Qc — 1, h = 1) and {gc — 5, 
/i = 3 and 5.5) for the solitons in the semi-infinite gap, and first and second finite bandgaps, respectively. 




(a) (b) 



FIG. 13: (Color online) Soliton-soliton collisions in the free space (e — 0) for Qc = 1, gd = —1, and fj, = —1. The 
initial velocities are c = ±0.2942 (a) and c = ±0.9651 (b). 



into a bound state, while, at high velocities, they pass through each other, as shown in Fig. [1^1 

The collision scenario is different in the opposite case, with the local attraction and nonlocal repulsion. 
As shown in Fig. [TH at small velocities the solitons bounce from each other. This feature is easily explained 
by the fact that the long-range interaction between the solitons is repulsive. The rebound is changed by the 
merger at intermediate values of the velocities. Finally, fast solitons pass through each other. 

Dependences of critical values of collision velocities, which separate different outcomes of the collision, on 
the strength of the local interaction are displayed in Fig. [151 
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FIG. 15: (Color online) Critical values of the collision velocity versus the strength of the contact attraction or 
repulsion for gd = 1 (a) and gd — —1 (b). In both cases, jj, — —1 and e = 0. 



In the presence of the OL, the sohtons can be made mobile (by application of a kick to them) if the 
nonhnearity in the model is weak enough; otherwise, the respective Peierls-Nabarro is very high. In the 
wcak-nonlinearity regime, the mobihty of GSs in the present model is quite similar to that reported in 
the local model with the self-repulsive interactions js^, as well as in the discrete model including the DD 
interactions I19l. 



VI. CONCLUSION 



This work presents results of the systematic analysis of one-dimensional bright solitons supported by 
contact and dipole-dipole (DD) interactions of opposite signs in BEG. In the absence and in the presence 
of the optical lattice (OL), stable soliton families have been found for the cases of local attraction and 
DD repulsion or vice versa. In particular, free-space solitons can be supported by arbitrarily weak local 
attraction if the DD repulsion is fixed; in the opposite case, there is a maximum value of the strength of 
the local repulsion, beyond which solitons do not exist (which was explained in an analytical form). In the 
model including the OL, a notable finding is a region of stability of subfundamental solitons (SFSs) in the 
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second finite bandgap. It is noteworthy too that, as seen in Figs. [TUl and [71 the gap sohtons (GSs) exist, in 
the case of the attractive DD interaction, if the contact repulsion is strong enough, and, in the opposite case 
of the repulsive DD interaction, GSs exist if the contact attraction is not too strong. 

Collisions between bright solitons in the free space were considered too. The collision scenario is the 
usual one in the case of the local attraction (merger and quasi-elastic passage at small and large velocities, 
respectively), while in the opposite case, when the local interaction is repulsive, a region of rebound was 
additionally found at smallest values of the velocities, which is explained by the long-range repulsion between 
the solitons. 
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